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Quantum phase transitions are usually observed in ground states of correlated systems. Remark¬ 
ably, eigenstate phase transitions can also occur at finite energy density in disordered, isolated 
quantum systems. Such transitions fall outside the framework of statistical mechanics as they in¬ 
volve the breakdown of ergodicity. Here, we consider what general constraints can be imposed on 
the nature of eigenstate transitions due to the presence of disorder. We derive Harris-type bounds 
on the finite-size scaling exponents of the mean entanglement entropy and level statistics at the 
many-body localization phase transition using several different arguments. Our results are at odds 
with recent small-size numerics, for which we estimate the crossover scales beyond which the Harris 
bound must hold. 


Progress in the study of interacting localization has 
brought to light new eigenstate phase transitions in the 
highly excited states of disordered isolated quantum sys¬ 
tems m- Across an eigenstate phase transition, prop¬ 
erties of individual many-body eigenstates exhibit singu¬ 
larities. Such transitions may be invisible in the tradi¬ 
tional equilibrium statistical ensembles, especially when 
they involve an extreme breakdown of ergodicity as in 
the many-body localized (MBL) phase [TH55]. Many- 
body localized phases have recently been experimentally 
observed in cold atomic and trapped ion systems [34H36] . 
The general theory of eigenstate transitions is largely un¬ 
known as the usual framework of the equilibrium critical 
theory need not apply. Thus, it is useful to find general 
constraints to guide their study. 

The most dramatic putative example of an eigenstate 
phase transition is the delocalization transition from a 
localized phase into a thermal phase satisfying the eigen¬ 
state thermalization hypothesis (ETH) [571 - IM| . This sig¬ 
nals the restoration of local thermalization at the level of 
individual eigenstates. To date, both numerical and phe¬ 
nomenological studies of delocalization in one dimension 
have found this kind of transition [HigEoiiioHia]. 

A widely used order parameter for this transition is 
the disorder averaged entanglement entropy [S'] in eigen¬ 
states, which shifts from an area law in the localized 
phase to a volume law in the thermal |19j . Assuming 
that the transition is continuous and that the entangle¬ 
ment entropy density [s] of a small fraction of the system 
satisfies a finite-size scaling ansatz, 

[s]{L,S)^^s{L^/''S) ( 1 ) 

we show that a = 0 and v > 2jd where d is the spatial 
dimension. We thus prove a Harris bound for the mean 
entanglement entropy density at the MBL-ETH transi¬ 
tion. This complements Grover’s recent result[28] on the 
critical value of the scaling function. The elementary 


arguments we present also apply directly to the mean 
level statistics ratio [r], another commonly used diagnos¬ 
tic [13]. 

The original Harris criterion that v > 2jd concerned 
the stability of clean equilibrium transitions to the in¬ 
troduction of disorder [33]. Chayes, Chayes, Fisher and 
Spencer (CCFS) generalized this exponent bound to arbi¬ 
trary disorder-driven transitions by defining a finite-size 
length ^ps associated with the tails of the distribution of 
the order parameter [331 |3S] ■ When a traditional coarse- 
graining picture holds, central limiting behavior governs 
the distribution of the order parameter away from the 
transition. We then expect ^ps to coincide with the cor¬ 
relation length governing the mean. This is the reason 
the CCFS result is often interpreted as a proof of the 
Harris bound. 

This argument could be applied to the entanglement 
entropy of the localized phase in the vicinity of a delo¬ 
calization transition in d > 2. Crudely, the entanglement 
entropy of a sub-region comes from independent blocks 
of size ^loc straddling the boundary of the region, where 
^loc is a localization length. Thus, for sub-regions signif¬ 
icantly bigger than ^loc, the entanglement entropy shows 
central limiting behavior. Consequently, the finite size 
scaling of the mean must be controlled by ^ps and obey 
the Harris bound. 

While this argument is appealing, it relies on several 
assumptions. First, it does not apply to non-spatial or¬ 
der parameters such as the level statistics ratio. Second, 
it doesn’t apply to d = 1 where most studies of the tran¬ 
sition take place. Finally, it ignores the possibility of 
Griffiths type effects that could spoil the central-limiting 
behavior. For example, the mean susceptibility in the 
paramagnetic phase of the random transverse field Ising 
model functions as an order parameter for the transition 
into the Griffiths phase. However, it does not have cen¬ 
tral limiting behavior and violates the Harris bound. See 
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also examples in Ref. [37]. 

The mean scaling bounds we derive in this article are 
independent of these considerations. In particular, we 
can derive Harris bounds in any d for [s] and [r] with no 
reference to detailed probability distributions. 

In the following, we first derive an elementary lemma 
regarding the scaling of the means of bounded random 
variables. This immediately leads to an exponent bound 
V > 2/(d + 2a) for any bounded order parameter with 
a scaling form as in Eq. Q. For completeness, we re¬ 
cap the result of CCFS and the associated dehnition of 
We then catalogue the behavior of several common 
order parameters used in the study of eigenstate phase 
transitions and discuss the relevance of these bounds to 
current numerical studies. We conclude with some open 
questions for future study. 

Lemma — Consider a collection of L‘^ sites subject to 
independent random disorder fields (or couplings) h with 
probability p\(h) dependent on disorder parameter A. 
For any bounded random variable Y S [—1,1], 


m 

d\ 




( 2 ) 


where = f dh —5— ( jg an 0{LP) constant 

px{h) \ d\ J 

and [•] denotes averaging with respect to the disorder. 

Proof — The lemma is a direct generalization of the 
proofs presented in |45L I46j . For simplicity, we review 
the argument for the case of a discrete bimodal distri¬ 
bution p\{h) = Xd{h — hi) -I- (1 — X)6{h — / 12 ). In this 
case, the probability of a disorder realization may be sim¬ 
ply expressed PrjH] = A""^(I — A)”^ where is the 

number of sites with held value hi in realization f2. We 
differentiate 

[r]=^Pr[H]y[H] (3) 

n 

with respect to A to obtain, 

(4, 


with a = yj A“^ + (1 ~ -^)~^ as required. The generaliza¬ 
tion to discrete distributions with q > 2 disorder values 
is straightforward while that to continuum distributions 
{q —>■ 00 ) requires somewhat more care [l^. QED. 

Taking Y to be an order parameter X, whose mean 
diagnoses a transition at d = A — Ac = 0, we immediately 
obtain the following theorem. 

Mean Theorem — Any bounded random variable X 
which satishes the hnite size scaling ansatz [35] . 

[X]{L,S)^^XiL^/’'S) (7) 


near a critical point at 5 = 0 has a scaling exponent 


2 

jy > - - 

d 2 q . 


( 8 ) 


This bound is not as strong as the usual Harris bound 
1 / > 2/d in the presence of a non-zero scaling dimension 
a. However, it assumes only the commonly used hnite 
size scaling hypothesis, while the following theorem, due 
to CCFS [35], requires a different hypothesis: 

Tail Theorem (CCFS)— Consider any random vari¬ 
able X which satishes the following two conditions with 
hxed 0(1) constants a and c > 0: 

• At the critical point at 3 = 0: Pr[A > a] > 2c for 
all L large enough. 


• Away from the critical point at d > 0: there exists a 
maximum length f,FsiS), such that Pr[A > a] < c 
for all L > fps- 

Then, the hnite-size scaling length is bounded below, 

fps > (9) 


where (3 = is an 0(1) smooth function of <5. 

Applications — Using the above results, we constrain 
the behavior of commonly used order parameters for sev¬ 
eral different eigenstate phase transitions. 

Entanglement Entropy — The entanglement entropy S 
of a subregion A of linear dimension Lji within an eigen¬ 
state \En) is dehned by 


S{L,LA) = -TrpA\ogpA ( 10 ) 


Taking the absolute value, recalling that |y| < 1 and 
applying the Cauchy-Schwarz inequality, 


dm 

dX 



(5) 


The radicand is simply a joint moment of ni and n 2 = 
L'^ — ni. Thus, 


dm 

dX 




( 6 ) 


where pA = Tr^ \En){En\. In the thermodynamic limit 
L —>■ 00 , the entanglement entropy density s = SjE^ 
functions as an order parameter for the delocalization 
transition: it scales to zero in the localized phase and a 
finite value in the delocalized as La —>■ oo. 

In disordered systems, the spectrum itself is random 
and so one needs a prescription for choosing the eigen¬ 
states of interest near an energy E. One choice is simply 
to take the eigenstate |A„) with the closest energy to E. 
Another possibility is to average over a subextensive en¬ 
ergy window around E. The results below regarding the 
entanglement entropy S are indifferent to this choice. 
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The entanglement entropy density s is bounded be¬ 
tween 0 and logq where q is the local Hilbert space di¬ 
mension. Assuming a finite-size scaling ansatz for [s] 
near a transition at disorder parameter 5 = 0, requires a 
two-parameter scaling form, 

( 11 ) 

where both and LaIL are held fixed as L —?► oo. 

The scaling ansatz allows for the exponents and scaling 
functions to depend on other system parameters [25] • For 
example, they could depend on the energy density e = 
E/L‘^, which we have also held fixed in this limit. 

At small La/L, we expect to access the scaling be¬ 
havior in the thermodynamic limit at La/L = 0. In 
the localized phase, the correction due to fixed La/L 
should be exponentially small in L as the contributions 
to entanglement come from the boundary itself. In the 
delocalized phase [S'] remains extensive although, in prin¬ 
ciple, the density [s] could be reduced at non-zero La/L 
from its thermodynamic value. This does not happen in 
phases satisfying ETH, as the correction to [s] (due to 
Page [SD]) vanishes as l/L'^. The argument below only 
requires [s] to be non-zero on the delocalized side and 
does not use this precise value. Thus from here on, we 
assume that La/L is held fixed to some small fraction as 
the thermodynamic limit is approached. 

There are two basic scenarios for the delocalization 
transition depending on whether the delocalized phase 
satisfies ETH. 

Localized-Delocalized (ETH) — In the ETH phase, the 
eigenstate ensemble reproduces the statistical ensemble 
for all local measurements [57II39| . In particular, the 
entanglement entropy density [s] is the thermal entropy 
density. This implies that [s] jumps at the transition 5 = 
0. The finite-size scaling hypothesis requires that a = 0 
as a scaling dimension a > 0 would impose a continuous 
turn on [s] ^ 6"°'. Thus, by the Mean Theorem, [s] must 
exhibit scaling with an exponent satisfying the Harris 
bound iz > 2/d. 

In d > 2, we argued before that the probability dis¬ 
tribution of s in the localized phase is likely Gaussian. 
If this is the case, the GOES length coincides with the 
finite-size scaling length controlling the mean. In d = 1, 
as less is known about the distribution of s, we merely 
speculate that the CCFS length can be constructed. Al¬ 
though it would also satisfy the Harris bound, it need not 
coincide with the finite-size scaling length controlling the 
mean. 

Localized-Delocalized (non-ETH) — If the delocalized 
phase fails to satisfy ETH, there are no further general 
constraints on the value of [sj. Consequently, the Mean 
Theorem provides a weaker bound as the presence of a 
non-trivial scaling dimension a is not ruled out. If the 
CCFS length exists, it clearly satisfies the stronger Harris 
bound. 


Unlike the MBL-ETH transition, there are several con¬ 
trolled examples of the localized-delocalized (non-ETH) 
transition. We review two as they illustrate the finite-size 
scaling behavior for [s]. However, both examples build on 
underlying critical theories with traditional coarse grain¬ 
ing descriptions and thus the mean scaling behavior is 
constrained by CCFS. 

The first example is the Anderson model wherein free 
fermions hop in a disordered background potential in 
d = 3[Tn]. As a function of decreasing disorder strength, 
this model undergoes a transition from a localized phase 
to a delocalized phase at infinite temperature. The delo¬ 
calized phase is non-ETH as a hnite density of the single¬ 
particle wavefunctions are localized at any non-zero dis¬ 
order strength. The value of [s] « p log 2 where p is the 
fraction of delocalized single particle states. 


The infinite temperature eigenstates are Slater deter¬ 
minants obtained by filling each single-particle orbital 
with probability one-half. Their entanglement structure 
follows directly from the spatial properties of the single¬ 
particle wavefunctions [STJ [52]. In particular, Eq. (11) 
should hold with the single-particle localization length ( 
as the diverging length scale 6~''. The d = 3 localiza¬ 
tion length exponent zz « 1.5 > 2/c? so that the Harris 
bound is satished IS31S1I- Note however that the scaling 
dimension a is non-zero at the transition because the ex¬ 
tended states are multifractal and have vanishing density. 
The scaling dimension can be estimated by assuming that 
the extended states contribute order one and the critical 
states contribnte ^ 1/^(4 to [sj. Using the density of 
extended states on the delocalized side p ~ 5^/"^ and the 
density of critical states ~ , we obtain 1/v -\-h = a 

and a = \/2v « 0.3. 


The second example is the random Clifford circuit 
model in d > 2 (see Ref. [55| for details). This is a 
tractable Floquet qubit model in which the localization- 
delocalization transition maps directly onto site percola¬ 
tion on the underlying lattice. Similar to the Anderson 
model, the delocalized phase possesses a finite density of 
fully localized conserved quantities which reduce the en¬ 
tanglement entropy density from log 2. In each disorder 
realization, S' in a region A is approximately determined 
by the number of sites within A which are connected 
to the boundary. As finite-size scaling holds near the 
percolation transition, Eq. 0 holds with the percola¬ 
tion length as the diverging length scale. The associated 
= Vperc satishes the Harris bound in all dimensions 
d>2 [51117]. The scaling dimension a is determined by 
the fractal dimension of the infinite cluster. 


Finally, we note that both examples demonstrate that 
the two-parameter scaling at fixed La/L accesses the 
same physics as that of the intrinsic limit La/L = 0. 


Level spacing ratio — The level spacing ratio r of three 
consecutive energy levels E„_i < En < En+i is defined 
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as: 


r„ = min 


f ae„ 


\AE, 


n+l 


AEn+1 \ 
AE„ ) 


( 12 ) 


where AEi = Ei — Ei-i. The prescriptions described in 
the previous section can be used to define r{E) near an 
energy of interest E. 

The disorder distribution of r distinguishes the ETH 
phase from localized phases [131158]. In ths ETH phase, 
the levels repel, exhibiting GOE statistics and [r] « 
0.527. In the localized phase, the levels exhibit Poisson 
statistics and [r] « 0.386. The value of [r] is an order pa¬ 
rameter for the MBL-ETH transition which jumps across 
the transition. The value at the transition itself is under 
contention; recent work [59j suggests a third intermediate 
value determined by semi-Poisson statistics. 

Clearly, r is a bounded random variable (0 < r < 1). 
Assuming that [r] obeys a finite-size scaling ansatz of the 
form in Eq. Q at the MBL-ETH transition, it immedi¬ 
ately follows that a = 0 as [r] jumps at the transition. 
Thus, u >2/dhy the Mean Theorem. 

The divergence of the length scale associated with 
the tails is more subtle as ^ps rnay not exist. In partic¬ 
ular, if the order parameter X is chosen to be r„ where 
En is the closest energy to E in each disorder realiza¬ 
tion, then its distribution has width of order one in ei¬ 
ther phase as L —>■ oo and ^ps is not defined. If, how¬ 
ever, we choose X to be the level spacing ratio averaged 
over a sub-extensive energy window AE in each disorder 
realization, then we expect the distribution narrows as 
L —>■ oo and S^ps can be defined. 

In the tractable examples of transitions between local¬ 
ized phases and delocalized non-ETH phases, [r] is not 
a good order parameter for the transition. In the non¬ 
interacting Anderson model, the many-body levels al¬ 
ways exhibit Poisson statistics and [r] « 0.386. Just to be 
clear, the level statistics of single-particle level spacings 
jumps at the transition. Similarly, in the Clifford model, 
the presence of strictly localized integrals of motion ev¬ 
erywhere in the phase diagram also leads to Poisson-like 
statistics. 

Relationship to current studies — In the last few years, 
studies of the MBL-delocalized transition have taken two 
approaches: numerical diagonalization of small systems 
and phenomenological real-space renormalization groups. 
The latter studies have found correlation length expo¬ 
nents consistent with the Harris bound, although their 
relationship to mean microscopic observables is unknown 

[miiiisD]. 

Recent numerical studies in one dimensional spin 
chains have reported scaling exponents violating the Har¬ 
ris bound jaiMi as measured by order parameters which 
satisfy the conditions of the Mean Theorem. This ei¬ 
ther signals the failure of the finite-size scaling hypoth¬ 
esis, or that the largest systems are still too small to be 
in the scaling regime. One can estimate a crude upper 


bound on the system size L* beyond which the observed 
scaling forms cannot hold. This follows as the Lemma 
provides an absolute lower bound on the width of the 
scaling window and not just on the exponents. Unfor¬ 
tunately, L* « 500 — 5000 for both numerical studies 
depending on which order parameter is considered. It is 
not surprising that L* is so large as 1) the Lemma uses no 
physical information about the order parameter, 2) the 
observed exponent likely increases with L, and 3) even in 
the scaling regime, Eq. ([^ is likely not saturated. Never¬ 
theless, L* provides some measure of asymptotic system 
sizes. We have estimated L* for the standard deviation 
CTs = \/[(s — [s])^], the mean [r], and the entanglement 
entropy generated in a local quench [AS]/L. All of these 
satisfy the conditions of the Lemma. 

We note that the constant a in the Lemma is formally 
infinite for the box distributions used in the numerical 
studies. In order to make the above estimates, we replace 
these distributions by Gaussians of the same width. It 
would be interesting to redo the numerical studies with 
disorder distribution which optimize the bound. 

Discussion — In this article, we have derived elemen¬ 
tary bounds on the finite-size scaling behavior of several 
order parameters associated with eigenstate phase tran¬ 
sitions. In particular, the Harris bound holds for [s] and 
[r] at the localized-ETH transition and a weaker bound 
holds at general delocalization transitions without fur¬ 
ther assumptions. Our results are complementary to the 
celebrated CCFS bound. 

Our mean results straightforwardly generalize to 
higher cumulants of bounded order parameters includ¬ 
ing the oft-studied variance of the entanglement entropy. 
They also can be readily applied to eigenstate phase tran¬ 
sitions within the localized phase between states possess¬ 
ing distinct long-range quantum orders 13 i- For ex¬ 
ample, the spin-glass order parameter and the spin-spin 
correlation function at half system size in the random 
transverse field Ising model both satisfy the mean the¬ 
orem. In this case, the stronger Harris bound can be 
proven using the strong-disorder renormalization group 
treatment [6T] . 

Finite-size scaling likely applies even if the delocaliza¬ 
tion transition is “first-order” [62] . At a first order tran¬ 
sition, there are no divergent length scales in the thermo¬ 
dynamic limit. Nonetheless, divergent finite-size scaling 
lengths control the scaling windows of [s] and [r]. These 
satisfy the Harris bound by the Mean Theorem. 

At the other extreme, there could be multiple diverg¬ 
ing length scales visible in the order parameters. This 
scenario would lead to the breakdown of the finite-size 
scaling hypothesis. It would be interesting to place gen¬ 
eral constraints on the transition in this case. 
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